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, In order to fix the parameters for predictions of hard photon and pion production 

|- in An + Au collisions at -y^ = 200 GeV, proton-proton and proton-nucleus data 

' are analyzed in perturbative QCD in the energy range -y^ fii 20 — 60 GeV. 
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1 Introduction 



Interest in direct photon production at RHIcB'i is motivated by the long 
^ I mean free path of photons. Hard photon production can be calculated in 

CO . the pQCD-improved parton model. Experiments need to separate the direct 

' photons from the tt" (and to^) two-photon background. It is estimated that 

the direct photon yield has to be at least 10% of the tt" yield to be measurable. 
I We have embarked on the following project: (i) fix the pQCD description of 7 

Q\ . and 7r° production in pp collisions; (ii) address nuclear effects in pA collisions; 

0^ I (iii) calculate direct 7 and 7r° production in AA collisions at ^/s — 200 GeV. 

fH , Our first results, which pertain to (i) and (ii), are summarized in Ref. 3. 

P^l Parton cross sections are calculable in fifiCD at high energy to leading 

order (LO) or next-to-leading-order (NL0).ai3 The parton distribution func- 
tions (PDFs) and fragmentation functions (FFs), however, need to be fitted 
. to data. In recent NLO calculations the various scales are optimized. □ Alter- 

natively, an additional non-perturbative parameter, the width of the intrinsic 
transverse momentum (fcy) distribution of the partons is introduced. Q We 
choose the latter method and use kr phenomenologically in pp collisions, ex- 
pecting that its importance will decrease in higher orders of pQCD. 



2 Proton-proton collisions 

In the lowest-order pQCD-improved parton model, direct pion production can 
be described in pp collisions by 

- E / ^^1^^ /i(^i.Q') f2{x2,Q') Aab ^ cd] , (1) 



abed' 
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where fi{x,Q'^) and f2{x,Q'^) are the PDFs of partons a and b, and cr is 
the LO cross section of the appropriate partonic subprocess. The K-factor 
accounts for higher order corrections. 13 Comparing LO and NLO calculations 
a nearly constant value, K « 2, is obtained as a good approximation of the 
higher order contributions in the pt region of interest. U In eq.(^ D{zc, Q^) 
is the FF of the pion, with Q = px/zc^ where Zc is the momentum fraction jaf 
the final hadron. We use NLO parameterizations of the PDEsllj and FFs 
with fixed scales. Direct 7 production is described similarly. □ 

To incorporate the intrinsic kx, each iritegral is extended to /cT-space, 



gikr) = exp(- 



dx d'kr gikr) f{x,Q'). 



We approximate gikx) as 



fc|,/(fc|,))/7r(fc|,). Here (fc|,) is the 2-dimensional width of the 



The 



kr distribution, related to the average kx of one parton as (k^) =t 4\{hir\ 

We applied this model to data from pp tt^X and pp 
calculations were corrected for the finite rapidity windows of the data. The 
Monte-Carlo integrals were carried out by the standard VEGAS routine. Il3 
We fitted the data minimizing ~ '^{Data — Theory)^ /Theory^ in the 
midpoints of the data. Fig. 1. shows the obtained fit values for {k"^) and 
a calculated {kr). The error bars display a = A^j„ ± 0.1 uncertainty in 
the fit procedure. The lines guide the eye and indicate, a linear increase to a 
common value of (fcr) ~ 3.5 GeV at ^/s ^ 1800 GeV.EZi 

These values of (fcr) provide a satisfactory description of the data up to 
energies of -^/s « 65 GeV. The data/theory ratios cluster around unity, with 
no systematic trend in the deviations (different experiments show different 



> 

o 



: 1 l„ 1 1 1 1 1 1 1: 
: p+p -> 7t +X " 


■- p+p -> Y +X 












1 1 III 





10 20 30 40 50 70 100 10 
(GeV) 



2.5 


1 1 1 1 1 1 1 1 1. 


2 












0.5 



+ 

1 1 1 1 1 


f, , , t 



20 30 40 50 70 100 
(GeV) 



Figure 1. Best fits for tfie widths of the intrinsic transverse momentum distributions in 7 
and tt" production data in pp collisions as a function of cm energy. See text for the lines. 
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slopes). B'tS The maximum deviations are on the order of « 50%. 

To estimate the px value where the direct photon production cross section 
surpasses 10% of the pion cross section in a pp collision at ^/s = 200 GeV, we 
use the \/s dependence of (fcr) from Fig. 1, and obtain the results displayed 
in Fig. 2. The threshold is reached in this approximation at pT ~ 6.5 GeV. 
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Figure 2. Extrapolated 7 (full line) and 7r''(x0.1) (dashed) production cross sections as 
functions of transverse momentum in pp collisions at ^/s = 200 GeV. 



3 Proton-nucleus collisions 

Having isolated the {k"^) already present in 7r° and 7 production in pp colli- 
sions at tlie present level of calculation, we turn to the nuclear enhancement 
observedll3 in pA collisions. In minimum-biased pA collisions the pQCD de- 
scription of the inclusive pion cross section is based on 



=E Z^'^^aW Jdxi,2d^kT,^, gi{kT,,b) 

abed 



92{kTM^l,Q')f2{x2,Q') K^Ri^!^. (2) 

dt TT^c 

Here b is the impact parameter and tA{b) is the nuclear thickness function 
normalized as J d^b tA(b) = A. For simplicity, we use a sharp sphere nucleus 
with = 2poy^B^^lP, where Ra 1.14^1/3 and pa = 0.16 fm^^. 

The standard physical explanation of the nuclear enhancement (Cronin 
effect) is that the proton traveling through the nucleus gains extra transverse 
momentum due to random soft collisions and the partons enter the final hard 
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process with this extra fcr.O We write the width of the transverse momentum 
distribution of the partons in the incoming proton as 

{k'^)pA = {kT)pp + C -hpAib) . (3) 

Here hpA{b) is the number of effective nucleon-nucleon collisions at impact 
parameter h imparting an average transverse momentum squared C . Naively 
all possible soft interactions arc included, but such a picture leads to a target- 
dependent Cn A more satisfactory description is obtained with a "saturated" 
hpA, where it takes at most one semi-hard [Q^ ^ 1 GeV^) collision for the 
incoming proton to loose coherence, resulting in an increase of the width of 
its kT distribution. This is approximated by a smoothed step function with a 
maximum value of unity. The saturated Cronin factor is denoted by C^""^ . 

As displayed in the left of Fig. 3, C"°* ^1.2 GeV^ gives a good fit 
of pA 7r°X data with A = Be, Ti and WH The right panel of Fig. 3 
illustrates how C*°* = 1.2 GeV'^ describes 7 production at ^/s « 30 — 40 GeV. 




Figure 3. Left: Cross section per nucleoli in pA — > tt'^X reactions {A = Be, Ti, W). We show 
Qsat _ ^ 2 GeV^ (full lines) and C"* = (dashed lines). Lower panel: data/theory on a 
linear scale for lib* pW collision. Right: Cross section per nucleon in the pBe ■yX reaction 
at two energies S3 Solid lines indicate C"""* = 1.2 GeV^, dashed lines mean C*"* = 0. Lower 
panel shows data/theory for = 31.6 GeV. 
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We interpret C^°* as the square of the typical transverse momentum im- 
parted in one semi-hard colhsion prior to the hard scattering. This picture and 
the energy dependence of C""^* need to be tested as functions of A at different 
energies, and in particular at y/s — 200 GeV, for RHIC AA predictions. 
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